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Nonexistence of classical diamagnetism and nonequilibrium 
fluctuation theorems for charged particles on a curved surface 
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Abstract. - We show that the classical Langevin dynamics for a charged particle on a closed 
curved surface in a time-independent magnetic field leads to the canonical distribution in the 
long time limit. Thus the Bohr-van Leeuwen theorem holds even for a finite system without any 
boundary and the average magnetic moment is zero. This is contrary to the recent claim by Kumar 
and Kumar (EPL, 86 (2009) 17001), obtained from numerical analysis of Langevin dynamics, that 
a classical charged particle on the surface of a sphere in the presence of a magnetic field has a 
nonzero average diamagnetic moment. We extend our analysis to a many-particle system on a 
curved surface and show that the nonequilibrium fluctuation theorems also hold in this geometry. 



Introduction - The Bohr-van Leeuwen (BvL) theorem 
states that the average magnetic moment of a classical sys- 
tem consisting of charged particles in the presence of an ex- 
ternal time-independent magnetic field is zero in thermal 
equilibrium [T]. The proof is in principle simple: Since the 
free energy, calculated from the canonical partition func- 
tion, is independent of the external magnetic field, the 
average magnetic moment, which is the derivative of the 
free energy with respect to the magnetic field, is therefore 
identically zero. 

At first glance the statement of the BvL theorem may 
appear disconcerting and counter-intuitive as, in the pres- 
ence of an external magnetic field, the charged particles 
undergo orbital motion. Thus, one could expect the sys- 
tem to have an average nonzero magnetic moment. But 
according to the BvL theorem, this is not the case. To 
understand the null classical magnetic moment physically, 
it is often pointed out that the boundary of a system plays 
a subtle role [TH!] . It is argued that the charged particles 
in the bulk undergo orbital motion which gives rise to a 
nonzero diamagnetic moment, but there is also a param- 
agnetic moment arising due to incomplete orbits of parti- 
cles which bounce off the boundary in a cuspidal manner. 
This paramagnetic contribution exactly cancels the dia- 
magnetic one so that the net magnetic moment is zero. 
Recently, based on this intuitive picture, Kumar and Ku- 
mar have considered a finite system in the presence of 
a magnetic field [2] where there is no boundary, such as 
a case of a particle moving on the surface of a sphere. 



These authors have claimed, from a numerical analysis of 
Langevin dynamics, that there exists a nonzero classical 
diamagnetic moment. It has been argued that the nonzero 
magnetic moment arises due to the avoided cancellation of 
the diamagnetic moment in the bulk and the paramagnetic 
moment having no contribution due to the absence of any 
boundary. 

However, in this paper, we demonstrate that the above 
argument is incorrect. To this end, we consider a parti- 
cle moving on a closed surface without any boundary in 
the presence of a constant magnetic field. We show an- 
alytically that, in the long time limit, the system which 
is governed by the classical Langevin dynamics is indeed 
described by the equilibrium canonical distribution. The 
case of a particle moving on a sphere, considered by Ku- 
mar and Kumar earlier [2J, is one special case of a class 
of systems we study here. Recently, Kaplan and Mahanti 
[5] have considered Langevin dynamics of a particle con- 
strained to move along a circle, but since in this case the 
motion of the particle is unaffected by the magnetic field, 
their analysis does not necessarily disprove the claim of 
Kumar and Kumar in full generality. However, motion on 
a sphere does depend on the magnetic field and therefore 
it is a first nontrivial case which will be studied here. Our 
results show that the previously claimed role of a bound- 
ary for having a zero diamagnetic moment is a misleading 
one because the average magnetic moment, as argued here, 
must vanish even for a finite boundary-less system which 
is shown to have a canonical distribution in thermal equi- 



p-1 



P. Pradhan and U. Seifert 



librium. In principle, one can understand the reason for 
vanishing of the magnetic moment as following. In ther- 
mal equilibrium, the probabilities for a particle to have 
velocity v and — v, at any given position r, are equal and 
therefore the average magnetic moment (q/2c)(f x v) of a 
particle of charge q is zero. 

We also generalize our results to a Newtonian many- 
particle system with charged particles moving on a curved 
surface of an arbitrary shape in the presence of a time- 
independent magnetic held. We show that Liouville's the- 
orem holds in the absence of any thermal noise and damp- 
ing. By adding those, the thermal distribution becomes 
the usual canonical one. At the end, we demonstrate with 
the Jarzynski equality [6] and the Crooks theorem [7 that 
paradigmatic cases of the nonequilibrium fluctuation the- 
orems (reviewed in [8TI10]) are valid even for a system with 
particles on a curved surface in the presence of an exter- 
nal time-dependent magnetic field and other nonconserva- 
tive forces. The fluctuation theorems in the presence of 
a time-dependent magnetic field in Euclidean space have 
been studied in [HUT!?] . 

Newtonian dynamics of a single particle - We start in 
the simplest setting where we first consider deterministic 
motion of a single particle of charge q and mass m mov- 
ing in the presence of a uniform time- independent external 
magnetic field B applied along the z-axis. Newton's equa- 
tion of motion can be simply written as, 



q 



mv = — (v x B) 



(1) 



where v the velocity of the particle and B — Bz, z be- 
ing the unit vector along the z-axis. In the presence of 
the constraint that the particle moves on the surface of a 
sphere with radius a, it is convenient to switch to spher- 
ical coordinate {r, 6} where the respective unit vectors 
are denoted as f, 6 and the constraint is expressed as 
r = constant = a. By equating force and acceleration 
terms along 9 and 6 directions, Eq. (jXJ) can be rewritten 
as 

ma{8 - 4> 2 sin 9 cos 6)0 = (— a<p sin cos 9)0 , (2) 

c 



maid) sin 9 + 206 cos 



-{—aOcosOU. (3) 

c 



We convert these two second-order differential equa- 
tions to four first-order differential equations in the 4- 
dimensional phase space {0, 6, vg, 
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using the relations 9 = vg/a and 6 = {v^ — 
a06 cos 0)/ (a sin 0) respectively in Eqs. ([2]) and ([3]), and 
writing the cyclotron frequency w c = qB/mc. 

The time evolution of the phase space density p, at a 
phase space point {0, 6, vg, v^}, is governed by Liouville's 
equation which can be obtained from the equation of con- 
tinuity in the phase space, 



dp 
dt 



d(pVe) , d(pV^) d( P V ve ) d(pV^) 
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The time-independent solution of Liouville's equation, 
dp s t/dt = 0, will be required later when we consider 
stochastic Langevin dynamics and deterministic many- 
particle Newtonian dynamics of a system in the pres- 
ence of a magnetic field. Since dVg/dO — dV<p/d6 = 
dV V0 /dvg — and dVy^/dv^ — — (cosOfasinO)vg, the 
time-independent solution p st is simply given by 
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where the numerical value of the microcanonical probabil- 
ity distribution p m i c follows from the normalization con- 
dition 



dA(6, 6, v e , v^pnicSimivj + v\)j2 - E) 
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(10) 



with J dA = a 2 J* sin OdO J^ n d6 dvg dv^ de- 
noting integration over volume in the phase space 
{9, 6, vg, v,p} and E the total kinetic energy of the par- 
ticle. Clearly the factor of sin (9 in Eq. (JSJ comes from 
the measure of an infinitesimal area on the surface of a 
sphere in spherical coordinate. The statement of Liou- 
ville's theorem that the phase space is incompressible un- 
der Newtonian time-evolution can now be expressed as 
dp/dt = where p = (p/s'mO). Note that the time- 
independent magnetic field does not perform any work 
(nor do the constraint-forces) on the particle, implying 
dE/dt — m(vgvg + v^v^) — 0, and using Eqs. ([6]) and (O, 
we get 

(vflV„,+«*V,J=0. (11) 

Langevin dynamics - Now we analyse the problem con- 
sidered in Ref. (2J, i.e., the classical Langevin dynamics of 
a charged particle constrained to move on the surface of 
a sphere in the presence of an external uniform magnetic 
field. The Langevin equations are obtained by adding a 
noise term and a friction force to Newton's equations, Eqs. 
© and ©, 



vg 

V4, 



7vg + y/ 2jk B T/mfg . 
7^0 + y / 2^k B T/mf^ > 



(12) 
(13) 



in addition to Eqs. (|4|) and ([5]), where 7 is the fric- 
tion coefficient, k B is the Boltzmann constant, T is the 
temperature, / is the Gaussian white-noise satisfying the 
fluctuation-dissipation theorem, (f a (t)f a '(t')) = S aa 'd{t — 
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t') and the subscript corresponds to the component along 
unit vector a or a! with a, a' = 9, 6. The probability 
density P{9, 6, vg, v$) is governed by the Fokker-Planck 
equation |13j . 

dp = d{v B p) _ d(y^p) _ gqyp) _ gqyp) 

dt 89 86 dvg 
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To solve the steady state of the above equation, i.e. 
the solution for which dP s t/dt = 0, we try the equilib- 
rium ansatz P s t{9,6,vg,v$) — const x p st x Pi(vg,v^,) 
where p st is the time- independent solution of Eq. © and 
Pi{v e ,v,p) = exp[-m(^ + v^)/2k B T]. Using dPi/dvg = 
— {■m/k B T)vgP\, dP\/dv,p = — (m/k B T)v^Pi and Eq. 
(fTTj) . the steady-state probability distribution for a par- 
ticle moving on the surface of a sphere turns out to be the 
equilibrium canonical probability distribution function, 



Pca.n(9,6,Vg,V r f > ) 
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where £ = m(v 2 + v 2 )/2 and the normalization condition 
being J dAp can = 1. The probability of finding the particle 
at any point of the sphere is uniform. According to the 
general theory of Fokker-Planck equations [13] , any initial 
distribution will relax to a unique stationary distribution, 
i.e., in this case the canonical one, in the long time limit. 
Thus one can see that the average magnetic moment is 
zero since (vg) — (v$) — with this distribution, and 
hence there is no classical diamagnetism. 

Langevin dynamics with an external potential - The 
above analysis can be easily extended to the case when 
there is an external potential U(r) at a position r due 
to conservative forces in addition to an external magnetic 
field. The Langevin equation is then written as 



mv — q(v x B)/c — VJ7 — jmv + \/2 r )mk B T f , (16) 

where / is the Gaussian white-noise satisfying the 
fluctuation-dissipation theorem. When the particle is con- 
strained to move on the surface of a sphere, Eq. (fH)]) can 
be rewritten by equating respective components of forces 
and accelerations along 9 and 6 directions, 



vg = V V( 



(VU)g/m - 1V g + ^2jk B T/mfg , (17) 
(W) /m - 1V<f> + ^/2jk B T/mf 4> , (18) 



where the gradient VU{r) is expressed in spherical co- 
ordinate with (W)<? = (l/a)(dU/d9) and (W) = 
(1 /a sin 6) (dU/d(p) respectively being the 9 and 6 com- 
ponents of the external conservative-forces. The full set 
of Langevin equations now consists of Eqs. Qj, ([S]), (fTT|) 
and (fT5)) . The corresponding Fokker-Planck equation can 
be written as 

dp _ d{VgP) djv^p) gqyp) gqyp) 

dt 89 d6 dvg dv^ 
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It is straightforward to show that the time-independent 
solution of the Fokker-Planck equation, dP s t/dt = 0, can 
be written as P s t = (sin0)p can where the canonical prob- 
ability distribution p can is given by 



p ca n{0,6,vg,v^) = (e 



_ ( e -S/k B T 



)/z 



(20) 



with the total energy £ = [m(v 2 + vfy/2 + U(9,6)] and 
the partition function Z = J dAexp(— £/k B T). Clearly, 
the probability distribution, and therefore the partition 
function, are both independent of the magnetic field. 

Newtonian dynamics of a many-particle system - Now 
we consider a classical Newtonian system consisting of N 
interacting charged particles, with mass m and charge q, 
moving on the surface of a sphere in the presence of a 
constant external magnetic field B(r) which may be space- 
dependent. Newton's equation of motion of the i-th par- 
ticle can be written as 

(mt/,). = (-Jyi x B{n))) & - (v Fi f/({n}))_ (21) 

where fi, Vi position and velocity of the z-the particle, 
(. . .)a denotes a component along the unit vector a with 
a = 9,6, and U({fi}) is the scalar potential containing 
both internal pair-potentials as well as an external poten- 
tial. Similarly as in the previous single-particle case of 
Eqs. (|4j)-<j7j) , Newton's equations of motion can be writ- 
ten in spherical coordinate in 4iV-dimensional phase space 
{9 i ,6 i ,vg t ,v c f li } as 



V 8i - v (i) 



v <f>i 



a sm f 



V 



(i) 



%=V«-i(V^({rl}) 
m 



**,=v«-i(^(W)) , 



(22) 
(23) 

(24) 

(25) 



where V 



(i) 



a[(cos 9i/ a 2 sm 9i)v 2 1 . + (lu c cos 9i/a)v t p i ] 
and (y n U({fi})) § = {l/a){dU{{9 l ,6 i })/d9 i ), 
Vv} = -a[{cos9 t /a 2 s\n.9 i )vg l v ( i li + {<jj c cos9 i /a)vg l ] 
and (v r - U({n})) . = (l/awn0i)(dU({6i,itn})/dttH). 

The time evolution of the phase space density 
p({9i,6i,vg i ,v ( p i }) in 4jV-dimensional phase space is 
given by the continuity equation, 



dp 
dt 
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Using dV^/ddi = dV^/dfa = dV&'/dvo, = and 

dVv^/dvfa — — (cosOi/asmd^vg^ it is straightforward 
to see that the time-independent solution of Eq. (f26|) . 
dp s t/dt = 0, is given by 



PatiiOufyuVe^Vfa}) =p mic JJsinfl,; 



(27) 



The microcanonical probability measure is 
Pmic = l/f2(J5) = const and n(i?) is the to- 
tal phase volume of a constant energy surface 
in 4./V-dimensional phase space i.e., fl(E) = 
H^j d^^ve^v^SiEHe^^ve^v^}) - E) 
where dA is the infinitesimal volume-measure as defined 
below Eq. (fTU|) and the total energy of the system is 
£ = \{m/2)YJ$ t + vl.) + U({0 Z , ((>.,})}. Also note that 
dp/dt — with p — p/(X\i sinflj), which is the statement 
of Liouville's theorem. Therefore, assuming ergodicity, 
an isolated system in equilibrium can be described by a 
microcanonical distribution where any time-averages can 
be in principle calculated from the uniform probability 
distribution on a constant energy surface. Since, in this 
case, two configurations with the same position coordi- 
nates and opposite velocities, i.e, {6i, (pi, V0i , and 
{9i, (pi, — vei, — v<j,i}, are equally probable, the ensemble 
(microcanonical) average of the total magnetic moment 

J2i(l/^ c )((^i x Vi)) 1S identically zero. 

Motion on an arbitrary curved surface - We now con- 
sider dynamics of a particle of mass m and charge q which 
is constrained to move on a closed surface of an arbitrary 
shape in the presence of an externally applied constant 
uniform magnetic field B = Bz. The position of the 
particle in a three-dimensional space is specified by a set 
of general orthogonal curvilinear coordinates {g M } where 
/i = 1,2,3 and the respective unit vector is denoted as 
q^. The constraint that the particle moves on a surface 
can be imposed, without loss of any generality, by taking 
q 1 = constant or equivalently q 1 — 0. Then a spatial point 
on the surface q 1 = constant is specified by the remaining 
two generalized coordinates, i.e., {q 2 ,q 3 }- An infinitesi- 
mal line element dl in the curvilinear coordinate can be 
written as 



{dlf 



3 

£ 

i 



g^dq^dq" , 



(28) 



where g^ is the metric tensor of the curved surface, rep- 
resented as a 3 x 3 matrix. Since the coordinate sys- 
tem is orthogonal, the metric tensor is diagonal, i.e., 
9iw — {h^Sfnv where is called the scale factor as- 
sociated with the curvilinear coordinate g M and 5^ v is 
the Kronecker-delta function. For example, in spherical 
coordinate considered previously, q 1 — r, q 2 = 9 and 
q 3 = <p and the respective scale factors are h r = 1, he = r 
and h<f, = rsin#. We also define a matrix element g^ v 
which is inverse of the 3x3 matrix representing the met- 
ric tensor g^, i.e., Y^s=i 9^ S 9Sv = ^- Since the met- 
ric tensor is diagonal, the inverse is also diagonal, i.e., 



g 1 " = (h fi )~ 2 6 fJ _v An infinitesimal volume element dV and 
surface element dA of the surface q 1 — constant can be 
written as dV = h\h2h^dq 1 dq 2 dq 3 and dA = h2h^dq 2 dq 3 . 

First we consider the deterministic motion of a single 
particle. Newton's equation of motion can be written in 
the covariant form If 4 , 



3 

£ 



El 

h x 



(29) 



where F x is the component of an external force along unit 
vector q x and the affine connection T x in an orthogonal 
curvilinear coordinate with g^ = S^g^ can be written 

as 



5,, 



dg 



xx 



dg. 



XX 



dq" 



iA 



(30) 



Note that we do not follow the Einstein summation con- 
vention here. Now defining the velocity v x along an 
unit vector q x as v x — h\q x and using the relations 
q x = (v x - h x q x )/h x , h\ = ^(dhx/dq^ in Eq. (22), 
we get a set of two first-order differential equations, 



(31) 



3 



dh x v x , A v» v v 

h x 2.^ T v 



fi—2 



hfj, h v 



— F — V a 

m A mag — v v x 5 



(32) 



where f^ag * s the component of the Lorentz force along 
the unit vector q x acting on the particle in the presence 
of a magnetic field. Note that the indices A, p and v 
are specified to take values 2 and 3 only as the motion 
along unit vector q 1 is not possible, i.e., v 1 = 0. The 
Lorentz force is given by F mag = {q/c)(v x B) where 
v x B — (v 2 q 2 + v 3 q 3 ) x Bz. Now taking the compo- 
nents of the Lorentz force tangential to the surface, we 
get F 2 lag = (qB /c)v 3 cos 9 and F^ ag = -{qB/c)v 2 cos 9 
where cos 9 = (q 1 ■ z), 9 being the angle between unit 
vectors q 1 and i. The angle 9 is a function of general- 
ized coordinates {q^}, i.e., 9 = Q({q^}). For example in 
spherical coordinate, using explicit expressions of the scale 
factor h\ and the affine connection T x , one could check 
that Newton's equations of motion in Eq. (|3T|) and ([32]) 
reproduce Eqs. (|4j)- dTJ) where in this case = 0. Since 
the external magnetic force and the constraint-forces do 
not perform work on the particle, using Eq. (|32l) . one can 
check that 



5>V = o = ]Tv^ 



(33) 



A=2 



A=2 



implying the total kinetic energy energy £ = 

;A=2( 



(m/2) ^ A=2 (w A ) 2 = const in this case. 



p-4 



Classical diamagnetism and the fluctuation theorems 



Alternatively, to derive Newton's equations of motion 
of a particle constrained to move on a curved surface, one 
can as well start with the Lagrangian L — K({q x , q x }) — 
U({q x }) where K = (m/2) J^ A (u A ) 2 is tne total kinetic en- 
ergy and U is the total potential energy where the equa- 
tion of motion is given by d(dL/dq x )/dt — dL/dq x = 
[13] . Using the afhne connection in Eq. (IBT)]) . one can 
check that Newton's equation in Eq. ([29]) is indeed same 
as the equation derived from the Lagrangian. 

The time evolution of the phase space density 
p({q x , v x }) in 4-dimensional phase space {q 2 , q 3 ,v 2 , v 3 } is 
governed by Liouville's equation, 



at ^ 



A=2 



d( P v q .) i d( P v v x) 

dq x 



dv x 



= 0. 



(34) 



A straightforward calculation shows that the time- 
independent solution of Eq. ([Mf is given by 



Pst(q 2 ,q 3 ,v 2 ,v 3 ) = p m i c h 2 h 3 



(35) 



where the microcanonical probability distribution p mic = 
1/Q,(E) = const and the phase space volume Sl(E) of a 
constant energy surface is given by Vt(E) = (J dA)S(£ — E) 
with dA = J1a=2 h\dq x dv x the infinitesimal phase space 
volume- measure and £ = (m/2) [(v 2 ) 2 + (v 3 ) 2 ] the total 
kinetic energy. Liouville's equation can be cast in a man- 
ifestly covariant form in terms of a scalar density variable 
p = p/yfg where g is determinant of the metric tensor 
[T6] . Then Liouville's theorem of incompressibility of the 
phase space can be shown to have a form dp/dt = 0. Note 
that Liouville's theorem holds even if one adds a position- 
dependent nonconservative force on the right hand side 
of Eq. (13"2"1) . since dv x /dv x = in this case. This will 
be required later to prove the nonequilibrium fluctuation 
theorems. 

The Langevin equation in the presence of additional ex- 
ternal conservative forces, derivable from a scalar potential 
U{{q x }), can be written as 



q X = V q > 



(36) 

v X = V.„a - ( W V /m - iv x + x /2 1 k B T/mf q , , (37) 

where (VU) q x = (l/h\){dU / dq x ) is the component of the 
conservative forces along the unit vector q x and f q \ is a 
white-noise satisfying (f q >~(t)f \>(t')) = 5\\<5(t — t'). The 
Fokker-Planck equation, as in the case of a single particle 
on a sphere in Eq. (|19[) . can be written as 



dP 
~dt 

A=2 



dv x 



A=2 

3 a 

T — 

\=2 



gqyp) | d{v v *p) 

dq x 
jv x P + 



dv x 
-yk B T dP 



dv x 



Pat = (h2hs) the time-independent solution of Eq. Q34p 
and the canonical probability distribution p C im({q X , v x }) , 



Pcan({gV A }) = (e 



_ ( p -£/k B T 



)/z 



(39) 



where £ — [(m/2)J2\(' 



,\\2 



U({q x })] and the partition 



function Z = U X (J h x dq x dv x ) exp(-£ /k B T). The av- 
erage magnetic moment can be shown to be zero since 
(v x ) — with this canonical distribution. Note that here 
we have considered under-damped Langevin dynamics on 
a curved surface. For over-damped Langevin dynamics on 
a curved surface and the corresponding covariant formu- 
lation of the Fokker-Planck equation, see [17j[T8] . 

In the many-particle case in a microcanonical set up, one 
can generalize Eq. ([27]) to a system consisting of particles 
moving on a surface of an arbitrary shape and obtain the 
microcanonical distribution, 



Pmic({qi,v x }) = l/Q(E) = constant. 



(40) 



The total phase space volume Q(E) of the constant en- 
ergy surface in 4iV-dimensional space is given by Q,(E) = 
Yli \ if (h\dq x )dv x }5(£ — E), where i is the particle index, 
£ = [(m/2)J2 hX (v x ) 2 + U({q x ,v x })} is the total energy of 
the many-particle system. Since the total energy £ is an 
even function of velocity v x , two configurations, one with 
{q x ,v x } and other with {q x ,—v x }, are equally probable, 
and this implies that the average total magnetic moment is 
zero for the microcanonical distribution. Adding dissipa- 
tion and a noise term satisfying the fluctuation-dissipation 
theorem to the deterministic Newton's equation of motion 
of the i-th particle, one obtains the Langevin equation 
which again can be shown leading to the canonical distri- 
bution in the long time limit. 

The Jarzynski equality and the Crooks theorem - We fi- 
nally discuss two remarkable relations called the Jarzynski 
equality (JE) [B] and the Crooks theorem (CT) [7], which 
involve fluctuations of work done on a system driven ar- 
bitrarily far away from equilibrium. 

First we consider the Langevin dynamics of a particle, in 
the presence of a time-dependent external magnetic field 
and nonconservative forces F nc which include a nonconser- 
vative electric field —dA/dt induced by the time- varying 
magnetic field B = (V x A), A being the vector potential. 
The Langevin equation (|36p still holds and, now adding 
the component of F nc along unit vector q x to Eq. ([37]) . 
we get 



v x 



v v x--(yu) q ,--r 

m 



1 

rn 



.A 

nc 



2-fk B T 



U ■ ( 41 ) 



(38) 



Proceeding as before, the time-independent solution of Eq. 
I|38p is given by the equilibrium ansatz P st = p s tp ca . n with 



Let us consider a process in a symmetric time interval 
—l~<t<T where F nc is nonzero only in a subinterval 
< t < t with r <C T so that, at t = ±T, the system 
is described by the canonical distribution as in Eq. ([39"]) . 
The work W done by the external forces and the heat 
Q transferred from the heat bath to the system can be 
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written as 



() I m (-71/ + v / 2jk B T/mf q ^ v x dt, (43) 



W = 



-T 

where a(t) is an external time-dependent parameter in 
the potential energy U which varies only in the inter- 
val < t < t. Under simultaneous reversal of time 
and the direction of the magnetic field, the total work 
done by the external forces is odd, i.e., W —> —W as 



t 



-t, v* 



and B 



-B (equivalently A 



—A) since F x c , which includes —dA/dt, does not change 
sign. The ratio of the probability Vf of a forward path 
{q x (t),v x (t);a(t),F x c (t)} to the probability Vr of a re- 
verse path {q x (— t), — v x (— t); a(— t), F x c (— t)} is given by 



Vf 
Vr 



exp 



= exp(-/3Q) (44) 



where J- x = [mv x 
have used Eqs. ((33 



(I3H and 



- F x 

__Jl 

m. 



of the system £ = [(m/2) 



,A\2 



|, p = l/k B T and we 
With the total energy 
U] , and by using Eqs. 
(J3S1) and (SP-dUD, we get AS = J{d£/dt)dt = W + Q 
where A£ the change in the total energy of the system. 
The probability distribution Pp(W) of work Wf done for 
the forward protocol {a(t) , F nc (t)} can be related to the 
probability distribution Pr(W) of work Wr done for the 
reverse protocol {a(—t), F nc (—t)} with the direction of B 
also reversed, 



P F {W)= p l V F 5{W F -W) 
forward paths 

-PQPi_ 
Pf 



£ p f Vn6(W R + W)> 
reverse paths 

= e^ w - AF )p R (W), (45) 



where the subscripts i, f denote the corresponding initial 
and final quantity respectively. The canonical probability 
distribution of the system p = exp[— j3(£ — F)] is given 
in Eq. (]3"9"]) with F = — (l//3)lnZ being the free energy, 
AF = Ff — Fi, and we have also used Eq. (l44l) . £f — £i = 
A£ = W F + Q and W F = -Wr. The relation in Eq. 
([45]) is the statement of the CT which has been shown 
here to be valid even for a particle on a curved surface. 
The JE, (exp(-/9W)) = exp(-/3AF), can be obtained by 
integrating the CT [7] . It is straightforward to extend this 
analysis to a many-particle system governed by Langevin 
dynamics on a curved surface. 

In a microcanonical set up, the JE and the CT have been 
recently proved for a system consisting of particles mov- 
ing in Euclidean space in the presence of a time-dependent 
magnetic field and other nonconservative forces. For par- 
ticles moving on a curved surface, the proof essentially 



follows from the fact that Liouville's theorem still hold for 
the system and the heat bath combined in the presence of 
a magnetic field and other nonconservative forces [12] ■ 

Summary - We have shown that the classical Langevin 
dynamics for charged particles on a closed curved surface 
in a time-independent magnetic field leads to the canoni- 
cal distribution in the long time limit. Thus the Bohr- van 
Leeuwen theorem holds even for a finite system without 
any boundary and the average magnetic moment for such 
classical systems is zero. Our analytical results disproves 
the recent claim by Kumar and Kumar 2 a that a classi- 
cal charged particle on the surface of a sphere governed 
by Langevin dynamics in the presence of a magnetic field 
has a nonzero average diamagnetic moment. We also show 
that nonequilibrium fluctuation theorems hold for a sys- 
tem consisting of particles on a curved surface in the pres- 
ence of a time-dependent magnetic field and other non- 
conservative forces. 
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